Abstract The new spinor-unit field representation of the electromagnctism 
[21j (with quark and lepton sources) is integrated via minimal coupling with 
standard Einstein gravitation, to formulate a Lagrangian model of the very 
early universe. A completely new solution to the coupled Einstein-Maxwell 
equations, with sources, is derived. These equations are generalized somewhat, 
but not in a way that violates any physical principles. 

The solution of the coupled Euler-Lagrange field equations yields a scale 
factor a(t) (comoving coordinates) that initially exponentially increases TV e- 
folds from a(0) ~ to a\ = a(0)e N (N — 60 is illustrated), then exponentially 
decreases, then exponentially increases to ai, and so on almost periodically 
(Oscillatory cosmological models are not knew, and have been derived from 
string theory and loop quantum gravity.) It is not known if the scale factor 
escapes this periodic trap. 

This model is noteworthy in several respects: {1} All fundamental fields 
other than gravity are realized by spinor fields. {2} A plausible connection 
between the unit field u and the generalization of the photon wave function 
with a form of Dark Energy is described, and a simple natural scenario is 
outlined that allocates a fraction of the total energy of the Universe to this 
form of Dark Energy. {3} A solution of an analog of the pure Einstein-Maxwell 
equations is found using an approach that is in marked contrast with the 
method followed to obtain a solution of the well known Friedmann model of a 
radiation-dominated universe. 
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1 Introduction 

1.1 Description of inflationary cosmology 

Inflationary cosmology [T01I141 H] is a widely accepted and largely successful 
theory of the physics of the very early universe. Recent observational data 
from important experiments such as WMAP [511711115] and the Sloan Digital 
Sky Survey [251151123]. and others support [51171127] the main predictions of 
inflation theory. Inflation employs a period of accelerated expansion at the 
beginning of the Universe to solve the monopole, horizon and flatness prob- 
lems. The simplest models of inflation are based on the slow-roll of a scalar 
inflaton field. The simple slow-roll inflaton potential model does not have the 
"exit from inflation" issue of the basic cosmological constant scenario, since 
the period of cosmic acceleration is followed by a radiation-dominated era, 
and then by a transient matter dominated phase [23] . After some mean time 
period, the inflaton excites Standard Model degrees of freedom and the hot 
big bang commences. Quantum fluctuations during the early time evolution of 
the inflaton field generate primordial curvature fluctuation spectra [17] that 
evolve into macroscopic perturbations whose time evolution leads to cosmo- 
logical structure formation. Most models predict a scale-invariant spectra of 
energy density perturbations and gravitational waves that act as seeds for 
structure formation. Moreover, the power spectrum of curvature fluctuations 
at horizon crossing during inflation may be related to the angular power spec- 
trum of cosmic microwave background (CMB) fluctuations at recombination, 
and yield predictions of the anisotropics in the cosmic microwave background 
radiation. 
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Yet, to date, no inflaton-type field model has gained universal acceptance. 
In this paper we study a new cosmological model based on the coupled {1} 
Einstein, {2} electromagnetic (with quark and lcpton sources) and the {3} 
unit field u equations that exhibits an oscillatory form of inflation/deflation. 
This model may also provide some insight into the nature of Dark-Energy. 
The photon field is represented according to [3T] using a real bosonic spinor 
field, and gives rise to the exact conventional interacting local [/(l)— invariant 
Maxwell equations after inflation has ended. "Maxwell" refers to the electro- 
magnetic sector of the Standard Model. For simplicity, we consider only one 
generation of leptons and quarks. Also for simplicity, this model employs a flat 
(k = 0) universe Friedmann-Lemaitre-Robertson- Walker (FLRW) line element 
with scale factor a = a(t) (comoving coordinates) for the Einstein sector. How- 
ever, we do not expect our qualitative results to change for non-zero, constant 
k, or with the inclusion of additional Standard Model degrees of freedom. 

The solution of the coupled field equations yields a scale factor a(t) = 
a[F(t)], where F(t) is defined in Eq.pQ below, that initially exponentially in- 
creases N e-folds from <x(0) ps to ai = a(Q)e N (N = 60 is illustrated; N is 
strictly a function of initial conditions), then exponentially decreases to a(0), 
then decreases to 0, then exponentially increases to ax, then exponentially de- 
creases to 0, and so on almost periodically. The energy density g, discussed 
below, determines the depth of the potential well of the periodic trap. It is not 
known if the scale factor escapes this periodic trap. The exponential change 
is with respect to F. The exact dependence is given by Eq.[2] and is only 
approximately exponential for large values of \F\. Nevertheless, for brevity we 
use the descriptor "exponential" . With respect to comoving time t the change 
is accelerated. The scale factor makes the transition from a(0) to a max in a 
finite comoving time that can be made as small as one chooses by changing 
initial conditions. Also, a max can be made as large as one wants and a(0) as 
small by changing initial conditions. 

Central to the new cosmological model is the unit field u [21], which is a 
fundamental spinor field realization of the unit element of the split octonion 
algebra. Throughout this paper we work with the components u a of u, 
rather than with u, to avoid the complexities of the non- associativity 
of the split octonion algebra. With respect to a canonical spinor basis of the 
algebra, its eight components u u a 1 a = 1,...,8, comprise a fundamental 
spinor that transform under a basic spinor representation of 50(4, 4). In this 
paper the unit field is a section u a = u a (x a ) a = 1, 4 of a certain spacetime 
fiber bundle whose technical specification is given in Appendix 1. The unit field 
may be employed to construct a local U(l)— invariant covariant derivative, 
enabling a basic (bosonic) spinor field representation of the electromagnetic 
sector of the Standard Model [3T] . 

In this paper the unit field is parameterized as 
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where F(t) = characterizes the Minkowski spacetime of Einstein's special 
relativity (but not uniquely). This Lagrangian model predicts that 

a{tf cosh [F(t)} = a(0) 3 cosh [F(0)] = constant, (2) 

i.e., a(t) 3 cosh [F(t)] is an exact integral of the motion. Since F(t) = charac- 
terizes the end of inflation, a( * enf ^ °£ (o y flat1oTl) = cosh [F(0)] 1/3 , where cosh [F(t)] = 

uu = [u 2 (t)} 2 + [u 6 (t)] 2 , is solely determined by the initial value F(0). (Here, 
the tilde denotes the transpose.) 

1.1.1 Radiation 

The Segre classification of the energy-momentum tensor [and, in virtue of the 
Einstein field equations, of the Ricci tensor] of a pure radiation field, which is a 
null Maxwell field, is the algebraic type A3[(ll,2)], with eigenvalue zero. There- 
fore the "radiation" that is modeled in the radiation-dominated Friedmann 
equations cannot be a solution of the Einstein-Maxwell equations because the 
Friedmann model stress-energy tensor is of type Al[(lll),l], which is that of 
a non-tilted massive perfect fluid with timclikc four-velocity. Type Al[(lll),l] 
is also different from the types Al[(lll,l)] (A cosmology), and Al[(ll)(l,l)] 
(non-null Maxwell field). In the Friedmann universe radiation is an incoherent 
superposition of waves with random phases and polarizations. If one integrates 
the stress-energy tensor I), „ of this electromagnetic radiation over solid angle 
then one finds that the average {T^u) is of type Al[(lll),l]. As an example, de- 
fine the lightlike Minkowski vector l a — (sm(#) cos^), sin(#) sin(<y9), cos(6>), 1) 
and use the definition of dT^ v given on page 521 of [55], which depends on the 
product Ftp. Then £ J Q 2 " ft l<* £ sin(9)d9dcp = diag (§, §, (T„„) 
defines/depends-on the radiation energy density g(t) and pressure p(t) func- 
tions, and appears in the Einstein equations G^ v = —2(T^^ that yield the 
radiation-dominated Friedmann equations. The equation of state relating the 
radiation pressure p and energy density g, after the averaging has been per- 
formed, iS P = g£. 

However, in this paper we solve an analog of the pure Einstein-Maxwell 
equations. We present a self-consistent model of gravity and radiation, with 
massive and electrically charged sources, that introduces the energy density 
g(t) and pressure p(t) |26j functions in order to make the total energy mo- 
mentum tensor homogeneous and isotropic (see Eq.|12|). There is no a priori 
equation of state relating p(t) and g(t) because p(t) is determined dynamically 
from the field equations (of course, g(t) is always determined from the field 
equations) . 

It should emphasized that we do not "second quantize" the quantum fields 
of massive, electrically charged fermions. A field equation for a massive and 
electrically charged fcrmion, information-equivalent to the usual Dirac equa- 
tion, is defined in terms of a real, eight component spinor field. The charge op- 
erator is not diagonal in this representation, as is well known from Schwingcr's 
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work with real eight component Fermion spinor fields [52] . But for simplicity, 
the real eight component Fermion spinor fields are not "second quantized." 
This simplifying assumption may be justified by the fact that we study a tem- 
perature regime far above the electron (resp., down quark) rest mass; therefore 
electrons and positrons (resp., quarks and anti-quarks) never annihilate effi- 
ciently. 

This paper is organized as follows: in the next section we define our notation 
and conventions; technical details are mostly relegated to the Appendices. In 
Section [3] the assumption that the universe is created in a non-equilibrium 
state is discussed, which is argued to imply that in order to make the total 
energy momentum tensor homogeneous and isotropic an additional term must 
be incorporated into the Lagrangian. Next the Lagrangian for this problem 
is defined. The Euler-Lagrange field equations are discussed in Section [5] A 
representative numerical approximation of the solutions of the field equations 
for the FLRW scale factor a(t), for the split octonion unit spinor field u a (t) 
and for the electromagnetic split octonion spinor fields are presented. Lastly, 
we summarize our results. 



2 Notation and Conventions 

Let G denote the Newtonian gravitational constant, c the speed of light in 

vacuum; is the Planck length I. We employ units in which AttG = 1 = 

h = c = eq = fj,Q. The elementary quantum of electric charge is e = v Ana ~ 
0.30282212..., where a is the fine-structure constant. The Planck time tp ss 
0.282095. For simplicity we assume that all Standard Model charges other 
than electric are zero, and neglect all contributions to energy, momentum and 
charge from electrically charged gauge bosons. Also, in this paper "leptons 
and quarks" means only electrons, and up and down quarks, and their anti- 
particles. 

Let (X3,i; g) be a spacetime. X3,i is a four-dimensional, Hausdorff, smooth 
manifold, and g is, almost everywhere [because the scale factor a(t) may peri- 
odically vanish] , a smooth pscudo-Ricmannian metric [a locally Lorentz metric 
of signature (+ H — h -) ] that is defined on X3.1. The manifold X3,i and the 
metric g are assumed smooth (C°°). We also require that X3 : i have a spin 
structure. We use the usual component notation in local charts. The tensor 
covariant derivative with respect to the symmetric connection associated with 
the metric g will be denoted by a semicolon. A partial derivative is sometimes 
denoted by a comma. 

In our model both the unit field and the photon field are realized as smooth 
sections of a spacetime spinor bundle. Technical details are discussed in the 
Appendices, and may be briefly summarized as follows: the unit field u, has 
8 real-valued field components u <H- u a = u a (x a ),a = 1,...8, a = 1,2,3,4, 
with respect to an oriented spinor basis of the split octonion algebra, that 
transform under the action of Spin(4,4,R) = 50(4, 4; R) as a basic type-1 8- 
component real-valued SO(4, 4;R) bosonic spinor field on X3.1. The photon 
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field <?2 <-> ^2" has 8 real- valued field components that transform as a basic 
type-2 8-componcnt real-valued Spin(4, 4, R) bosonic spinor field on ~&3 t i, 
ip2 = ^2 [ xOL )- [Upon restriction to SO(3, 4; R) the u a and transform in the 
same way.] Following Dirac, ip denotes the transpose of a (real) 8-component 
spinor -0. 

The new definition of the photon wave function is based on a new local 
f7(l)-invariant Lagrangian formulation of Maxwell's theory [2T| that is cast in 
terms of the unit field u and tp^ > an d which revises the definition of the pho- 
ton wave function, but whose associated Euler-Lagrange equations imply the 
conventional Maxwell equations exactly, in the flat spacetime limit. It cannot 
be overemphasized that in this formalism the electromagnetic interaction is 
carried by a real, 8-component, bosonic type-2 spinor field $2 O ip2 a , which 
defines the four-vector potential A a in the standard Maxwell theory via an 



interaction with the unit field, A Q 



ucr T a $2\ the t q matrices arc 



defined in the Appendix in Eq. [54] . In addition, the association to electric and 
magnetic fields (E, B) employs a Lagrange multiplier W\ O ip\ a [which is a 
real, 8-component, (bosonic) type-1 spinor field], related to the conventional 
(E, B) Minkowski frame components according to 
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(3) 



when u is given by Eq.[T]. (The field equations guarantee that 5( 4 ) = 
%)•) 



3 Non-equilibrium birth 

In this model we assume that at the onset of creation at t = the par- 
ticles/fields created by the first Big Bang are not in thermal equilibrium. 
We resolve the charge density of "all" electrically charged matter as platter = 

P%atter(t) = M JZo [(-iWectronft + HW™^ ^) + ( + f Kp(0 ^)] d 3 ^- 

Here p S p ec j es (i; v) = / ^ Pspecies^' x)d x denotes an associated en- 
ergy density distribution function and the integral is over the hypersurfacc 
t = 0. Let p® nti _ matter = Panti-matterit) denote the charge density of "all" 
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electrically charged anti-lcptons and anti-quarks (in this paper this means 
only positrons, and anti-up and anti-down quarks ). We assume that platter + 
Panti-matter = 0, so that the newly created universe has zero net charge (in 
fact, one of the Euler-Lagrange equations demands that this net charge be 
zero). However we assume that while the spatial components of the sum of the 
linear momenta of all particles and fields sum to zero, at comoving time t = 
a non-vanishing average current drift velocity J d instantaneously exists that 
satisfies 



-tie 



1 2 

(-l)p e lectron(t; If) + (- -) Pdown{t\ v) + (+g)Pup(*5 v) 



1 2 

(-typpositronfav) + {--)Panti-down{t;v) + (+-)panti-up(t',~v) 



d 3 lt. 

(4) 



It is assumed that the comoving z = x 3 -axis is aligned with the average 
drift velocity at t = 0. (When the field equations are solved it is found that this 
model predicts that very shortly after J = 0.) 7^ ^ sources a 

non-isotropic contribution to the total energy momentum tensor near t = 0. In 
order to make the total energy momentum tensor homogeneous and isotropic 
this non-isotropic term must be canceled by an additional non-isotropic term 
(see Eq. ) that must be incorporated into the Lagrangian. 

In passing we note that in the 1930s, Robertson and Walker independently 
proved that there are only three possible spacetime metrics for a universe 
that is homogeneous and isotropic. A spacetime having one of these three 
spacetime metrics must be associated with an energy-momentum tensor that is 
homogeneous and isotropic. Coley and Tupper [5] have shown that a fluid with 
an axial velocity may possess a homogeneous and isotropic energy-momentum 
tensor, and therefore source a Robertson and Walker spacetime geometry. In 
this paper the total energy-momentum tensor is homogeneous and isotropic, 
although its constituents are generally not. 



4 Lagrangian formulation of the problem 

Let E = det (e^J = \J—g 7^ 0, R denote the Ricci scalar and £Ei ns tcin-Hilbcrt 
KR. We write the total Lagrangian density for this problem as 

^ = 2 8tt G ^ ^Emstein-Hilbcrt + £u) 

+ £-EM + £-Fermions + £gp- (5) 

Here £ u is the Lagrangian for the unit spinor field u a , C em is the Lagrangian 
for the electromagnetic field ^Fermions is comprised of six Lagrangians 
that govern the dynamics of the massive and electrically charged electron, 



s 



Patrick L. Nash 



up quark and down quark, and their antiparticle fields. These holds are not 
second-quantized. C gp is somewhat similar to the familiar Lagrangian for an 
effective fluid that is characterized by an effective proper energy density g(t) 
and an effective pressure p(i). In this model both the energy density g(t) and 
the pressure p(i) are determined dynamically by the field equations. 



4.1 Lagrangian for the unit field u a 

Let 7° = E"^ T 1 -^ t^K The Lagrangian for the unit field u a i 



is 



£ U = E^= 1 5cTr (5 V^7 a u, a . (6) 

2V47tGuo-u 1 

The the a and r matrices are defined in the Appendix; a representation of the 
8 x 8 a matrix is given in Eq. [50] and a representation of the r matrices is given 
in Eq. |89j . The vertical bar denotes spinor covariant derivative (see Appendix 
E]). We note that t he terms t( 5 M 6 ) (and t^ 7 \ below) break the 50(4,4; M) 
symmetry down to 5*0(3, 1;K) . 



4.2 Lagrangian for the electromagnetic field ^ 

The interaction with the electromagnetic field possesses local U(l) invariance 
and incorporates the covariant gradient operator 



D a = I 8 x 8 V Q - t c ( ) -^=gt c . (7) 

\\Jwa w J \ a vuau 

The local [/(l)-invariant electromagnetic Lagrangian Cem for this problem is 

E2 



^ a P ± r« \ t c [ -== ) —=ot c I <P 2 

,VUCTU/ \ a Vucru J 



u 



(8) 



This Lagrangian gives the conventional interacting local [/(l)— invariant Maxwell 
equations in a fiat spacetime |21j . 
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4.3 Fermion $ Lagrangian 

For 3? = {fields associated to the electron, up quark and down quark, and 
their antiparticles}, there are terms in the total Lagrangian of the form 



£/•'< , »,„.,. [*] = -E $ cr 



Ir (5) T (6) T (7) # 



3>ot M K ""•^'fr-l r (5) T (6) r (7) * 



VtKTU 



(9) 



Here £ = characterizes the electric charge carried by the fields, and 

takes values = — 1, |,— | and | for <I> = (e,u,d,e,u,d). The charge 
operator is not diagonal and is given by 



g = er (5) T (6) T (7) 

The mass parameter M takes values M = M. e ,. 
and M-t = M^. The Fermion Lagrangian is 



(10) 



Ferrnicm 



A m (C F [*e] + U [*u] + CF [*d]) 



(ii) 



Here A m = A matter and yl a = A anti - m atter are order parameters that keep 
track of matter/anti-matter contributions to the field equations. 

In a flat Minkowski spacetime the Euler-Lagrange equations associated to 
this Lagrangian imply the same mass/four- momentum dispersion relationship, 
and minimal electromagnetic coupling, as does the conventional Dirac equa- 
tion. 



4.4 Density and pressure terms 

Our approach to implementing the consequences of the assumption of the non- 
equilibrium birth of the universe discussed in Section [3] , with its preferred (but 
short lived) initial direction, is to employ 



where 



/00 \ 

0p(i) 
\0 g{t) ) 



rpOL 



(12) 
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in the Lagrangian. Here g(t) and p(t) denote the energy density and pres- 
sure parameters, but they are not identified with a fluid. There is no a pri- 
ori equation of state relating p(t) and g(t) because p(t) is determined dy- 
namically from the field equations. g(t) may be computed as follows. We 
solve the field equations for a{t) [= DA(a(t),F(t), g(t), . . .), say] and a(t) [= 
DDA(a(t), F(t), g(t), . . .), say] and require that f t DA(a(t),F(t),g(t),...) = 
DDA(a(t), F(t), g(t), . . .) when the field equations are satisfied. g(t) appears 
linearly in the last equation, and can always be solved for. 

It should be emphasized that the total energy-momentum tensor for this 
problem is homogeneous and isotropic, although the constituents such as 
of Eq.[T3] are generally not. 

A possible alternative approach to employing Eq. [T^] in the total Lagrangian 
C is to include a self-interacting scalar field T in C using 

C r = -E I \g aSi d a T dpT + V(T) j . (14) 

Let St denote the action of the scalar field. The energy-momentum tensor for 
the scalar field is 

T<T> = -^=/^7 = d,Td v Y -g, v ^-g^ da Td t3 T + F(T)J . (15) 

Assuming the FLRW metric for g a p and restricting to the case of a homogc- 
ncous field T(t), this energy- momentum tensor T^J takes the form of a perfect 
fluid, with the energy density and pressure given by 

Q=\f 2 + V(T) 

V=\t 2 - V(T). (16) 
Therefore, unfortunately, this method cannot lead to a cancelation of a non- 

(T) 

isotropic term in the total energy-momentum tensor, because T^J is homoge- 
neous and isotropic. 



5 Field Equations 

All of the variations of the Lagrangian with respect to the fields arc first per- 
formed for the most general forms of the field functions, and then the special 
functional forms u a = ^(0,6-^,0,0,0,6^,0,0) , = diag(a(i), a(i), a(t), 1), 
and so on arc substituted into the resulting Euler-Lagrange equations. Unfor- 
tunately the full set of field equations is too lengthy to print (please see below 
for an exception) . However one combination of variational derivatives is simple 
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and important: 



5C | c F(t) 6C 



8u 2 {t) ' Su 6 (t) 

a{t) 2 (3 cosh [F(t)]a(t) + a(t) wah[F(t)] F(t) 

1 Fin d 



: e 



{a(t) 3 cosh [F(t)}} 



8V% dt 

= (17) 

Summarizing results, this model predicts that 

a(t) 3 cosh [F(t)] = a(0) 3 cosh [F(0)] = constant. (18) 

F(t) = is one characterization of the Minkowski spacetime of Einstein's 
special relativity. Since F(t) — > characterizes the end of inflation we conclude 

that a( * end °I(j," flat ' ion) = cosh [^(0)] 1/3 - F(0) ks ±3 x TV implies that inflation 
lasts for about N e- folds. 

The Euler-Lagrange equations imply that the magnetic field B vanishes 
but the electric field E has non-vanishing z-component (that exponentially 
vanishes) . 



5.1 Special case 

For the special case without Standard Model degrees of freedom, obtained by 
setting to zero the Fermion degrees of freedom and the electromagnetic ^ and 
<£i spinors, the Einstein sector of the field equations is 



SEX 



(/') 



diag 



/ a(t) [a(t){a(t)p(t) + g(t)} - a"(t)} ±a'(i) 2 \ 
a(t) [a(t){a(t)p(t) + g(t)} - a"{t)\ - \a'{t) 2 
a(t) [a(t){2a(t)p(t) + g(t)} - a"{t)] - \a'{tf 
p(t)a(t) A + 2g{t)a{tf 



%a'{tfa{t) 



With 



g(t) = sech[F(t)}, 

the field equations imply that 

P (t) = 
a(t) = 

o(t) = 

F(t) = T^VSgJt) coth(F(i)) 
&g(t) = ±2y/3g(t)g(t), 



±-a(t) VM*) 



(19) 

(20) 

(21) 
(22) 

(23) 

(24) 
(25) 
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where Eq. [T7] and Eq. [55] have been used to arrive at Eq. [5J] . Differentiating 
Eq. [55] with respect to t and setting the result equal to Eq. [53] yields 

g(t) = T^g(t) 3 ' 2 . (26) 

The field equations Eqs.[55],[55] and [5B] imply that 

g(t) = 9 Q (t)g(t) 



Q(t) = ^Q{t? 



(27) 
(28) 



5.1.1 Integrals of the motion 



Dividing Eqs.[55J and [55] by F(t) yields 



JL a (F) = ~a(F) tanh(F) 
JL e {F)= 1 - e {F) tanh(F) 



(29) 
(30) 

(31) 
(32) 

where Fi, ai and Qi denote initial values. These relations imply that g 2 a 3 = 
constant. 



Integrating, two independent integrals of the motion obtained: 
a 3 cosh (F) = constant = a 3 cosh (Fi) 
g 2 sech(F) = constant = g 2 sech(i 7 i), 



5.1.2 Integration of the field equations Eqs. 



W and 



Eq.[H] and Eq.[35] imply that g(t) = ±2g(t) 3 / A y/Tg i y / sech (Fi). The dif- 
ferential equation for F(t) has a regular singular point at F(t) = and is 
not analytic at F(t) = oo. Consequently g(t) is discontinuous at times when 
g(t) = 1, and switches Riemann surfaces when g(t) crosses zero. 

Shifting the comoving time origin from zero, the equations for a, g and g(t) 
may be easily integrated and the solutions expressed as 



a = Qi y sec\i(F) cosh(Fi) = ai 

g = gi Vcosh(F)sech(F, ; ) = g t , / — 

V 9 



(33) 
(34) 



git) 

sech(F(t)) 



scch (Fi) [t^- - ^cosh (Fi 



scchfF,) t 



{/ cosh (Fi 



= g t + 4 



/cosh(Fj) 
9 ft 2 



< t^SL < ^/ C osh(F l ) 
-^/cosh( J F 4 ) < t^p^ < 

(35) 
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where F{0) = ^ g{0) = 1 and Fi ( ^ F{0) in general), ft, di and Qi denote 
initial values. Therefore, without Standard Model degrees of freedom, F{t) is 
periodic in t with period 



I cosh. {Ft) 



9ft 



and angular frequency 



9ft 



cosh(F !; ) 2 V 



ft 2 ft 



(36) 



(37) 



that depend on initial conditions. The integral of the motion g 2 sech(F) = 
g 2 g = constant = § (^jr) ■ Note that the scale factor will make the transition 



from a,i to a r , 
than ^T. a,„ 



a(0) 



a, ^/cosh(i 7 i) = in a finite time that is less 

can be made as large as one likes by reducing ft, and T can 
be made as small as one chooses by increasing ft 2 ft. To address the horizon 
problem in a model that includes material, all that is needed is for a to be 
small enough for enough time for the system to come to thermal equilibrium. 
The fact that a vanishes guarantees this. Even if we can find a way to keep 
a nonzero, the fact that a oscillates over many periods will guarantee that 
ultimately a system is brought to thermal equilibrium. 
Alternatively, to solve the field equations we may define 

t 



r{t) = l vmm) 

a{t) = A[r{t)]G[r{t)} 4/3 
g{t) = G[r{t)} 4 
Q {t)=R[T{t)]G[T{t)]- 2 . 



(38) 



The field equations Eqs.[22],E5] an d 
single Riemann surface, are locally 



], re-expressed as functions of r on a 



^ GW = ±2 ,^( T ) = = A R(T) . 



(39) 



This special case is unphysical because of Eq.[3], J d(0) 7^ if. Nevertheless 
this case demonstrates that the physical F{t) may be almost periodic in t. 



5.1.3 Lagrangian for the special case 

A Lagrangian Lo may be defined for this special case. Let 
1 



L„ = 



a{t) - I a{t) V3g{t)^j (a{t) + | a{t) y/ 3 g{t) 



+ (g{t) - 2g{t) V 3 (?(*)) (g{t) + 2g{t) ^Jgjt) 



14 



Patrick L. Nash 



(40) 



Clearly Lo = when the field equations Eqs.[2i],[2S] and [55] are satisfied, 
which is the unique absolute minimum for this functional. The Euler-Lagrange 
equations associated to Lo are satisfied when the field equations Eqs.[2"2"],[25] 
and [2S] (and their higher derivatives, Eqs.[27] and [55]) are satisfied. Put x 3 = 
(a,g,g), j = 1,2,3. The canonical momenta are pj = ^|jLo. The Jacobian 

dSdxi * s non-zero, so we can solve for the x 3 as functions x 3 — v 3 (t,x 3 ,Pj) 
of the canonical variables (t,x 3 ,Pj). 

The canonical Hamiltonian function is given by M(x 3 ,pj) = — Lo + pjx 3 , 
and in this case is independent of t. Here it is assumed that the x 3 , on the 
right-hand-side of this definition are the inverse functions x 3 = v 3 (t, x 3 ,Pj). 
Lo is of class C°° if we exclude the hyperplane g = in (a, g, g,p a ,Pg,Pg) — 
phase space. It follows that the pj are continuously differentiablc functions of 
{t,x 3 ,x 3 ) in the same region. However the inverse functions v 3 (t,x 3 ,pj) and 
H(i, x 3 7 pj) are not continuous everywhere in phase space. These functions 
have branch points and cuts in the six-dimensional phase space, in virtue of 
the roots of the combined powers of the momenta and coordinates that arise 
in solving for the x 3 . For example, when a > and < p a << 1 then 

(8+8i\/3)a 2 -2(-l) 2/3 2 1/3 (-27p a +V729p a 2 -256aV) 2/3 q„ , , , 

a = v ,, 3 '- — w f^r- + \aJ3g] also 

6 2 2 /3 (-27p a + v /729p a 2 -256a 6 e 3 J Ha g i 

g m 24 ^ 2 - + 2 g\/3g, g ~ — \/3gg, which are real in the appropriate region 
of phase space. These relations are reminiscent of a charged particle moving 
in an external vector potential proportional to v / 3g(|a, 2 g, — g) with respect 
to some curvilinear coordinate system. Near the origin of momentum space 

M 3 Po 2 P\ P g 2 (2 \ 

'Zq ~aPa + 2gp g -p e g). (41) 



u \Qa?g 12p 3 48g 2 g v e \ 3 

Once the Hamiltonian Ho has been defined and operator ordering problems 
have been resolved one may in principle quantize the system and study the 
cigenstates of Ho. An understanding of the quantized system may aid in the 
interpretation of the periodic solutions of the field equations. 



5.2 Full problem 

For the full problem we find that F(t) has a regular singular points at F(t) = 
and is not analytic as F(t) — > oo. After one fixes the initial values at t = of 
the Fermions and the electromagnetic ^ and <Zq fields, with g(0) > 0, then 
our numerical experiments verify that g(t) > and that cos ^[jr(t) ] ~ C osh [ F (o) ] ' 
The rate of change of g(t) when a > compared to this rate when a < 
depends on A matter and A anU - m atter, as follows (using Eq.JTJj]): 
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Q{t)\a>o = -Q(t)\a<o + V2ipl{t)e- 1 { 

1 . 



A, 



1 



An I -IVlrf 



'¥2t)r i QW d {t) 



Wjt) T 4 QW e {t) 



<Pu(t)r i Q<P u (t) 
C(t)r 4 QC(t) 



A, 



r 4 r 3 ^(t)J [* B (i) r 4 r 3 # e (t)J + - r 4 r 3 «2>„(i) 

2 r- 



<Z> rf (i)r 4 T 3 <Z> rf (t)J + [<2> e (f)r 4 T 3 # e (i)J - - [<2>„(i) r 4 r 3 * u (t) 



The physics governing the escape of F and a from the periodic trap has 
not been discovered. If the value of the energy density g(t) decreases to the 
point where 

g(t) 



lim 

F(t)->o F{ty 







then limpn^a F(t) = 0, as well as limp( t )_yo d Jt^ = 0, n = 2,3,..., and 
therefore F — > 0. a is controlled by F [because a 3 cosh (F) = constant] 
and will escape the periodic trap and very rapidly transition to a constant 

°(*end of inflation) ~ «1» sincc F W = v '* > ^end of inflation- However the 
governing physics may have more similarity with quantum mechanical tunnel- 
ing through a barrier than the decay of the controlling parameter g(t) to zero. 
The investigation is ongoing. 

In numerical experiments we seek a solution of the form u a = ^= e ^ " 

(with F(t) = characterizing the Minkowski spacetime of Einstein's spe- 
cial relativity, but not uniquely), V%(t) = {Sf + <Jf e~ F(t) } an d = 
{5% - 5g e F W} V'i(i)- The fermions are assumed to have vanishing components 
<Z> a for a = 3, 4, 7, 8. We find that 



j(t) = ±g(t) V6 



(42) 



2e(*)-V?(*)' 

^ (- 2 - (^(o) 2 + (*sw) 2 + (*sw) 2 - 3 (*i(*)) 2 - 3 (m) 2 

+ 3 (#»(*)) 2 + 3 (*»(*)) 2 + 2 ((**(*)) 2 + (^( t )) 2 - (#»(*)) 2 - (*•(«)) 2 )) 

+ A m ((*&)) 2 + (*g(t)) 2 - (*§(*)) 2 - (*g(t)) 2 + 3 (ffj(t)) 2 + 3 (<2> 2 (t)) 2 
- 3 2 - 3 2 - 2 (*i(t)) 2 - 2 (<2> 2 (*)) 2 + 2 (*£(*)) 2 + 2 (#«(t)) 2 )) 

+ yl ro (M d (#3(*)*§(t) - 

+ M e (^ e m 5 e (t) (*)*«(*)) +M U ($Km 5 u (t) (*)*£(*))) 

+ 4, (M d (^(*)^(*)-*i(*)*SW) 
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where ± according to a ^ 0. 

Initial values are o(0) = 10~ 24 , F(0) = -180, £-(0) = 10 40 ,^(0) = = 
V> 2 (0), <P e (0) = ^-(1,0,0,0,1,0,0,0) = # d (0) (and their antiparticles) and 

U (O) = v/2^(l, 0, 0, 0, 1, 0, 0, 0) = i^(0). Figures [4-8] summarize the results 
of representative numerical integrations of the coupled Held equations. We 
note that experiments with ^1(0) ~ 10 4 ° ~ {V'i(O)} have been performed 
using arbitrary precision floating point numbers to avoid over/under- flows near 
t = 0. Calculations employing 450 bits for each floating point number seem 
to maintain accuracy. However, to date, high precision calculations (running 
several days) have only evolved the system to t = 10~ 98 T, where T is the 
expected period and is of order 1. 



6 Conclusion 

We have studied a new cosmological model that is based on a Lagrangian 
and couples Einstein, electromagnetic (with quark and lepton sources) and 
unit u field [3T] sectors. We have focussed on a particular solution of cou- 
pled Euler-Lagrange field equations that employs a flat (k = 0) universe 
Fricdmann-Lcmaitre-Robertson- Walker (FLRW) line element with scale fac- 
tor a = a(t) [comoving coordinates] for the Einstein sector. The solution 
of the field equations yields a scale factor a(t) that oscillates almost peri- 
odically, and admits the first integral a(t) cosh[F(t)] = a(0) cosh[F(0)] = 
constant. Since F(t) = OVi > t en( ^ f mna tion characterizes the end of infla- 
tion a( *cnd of inflation } = cosh [^(O)] 173 , where cosh [F(t)\ = uu = [u 2 (t)] 2 + 

[u 6 (t)] 2 , is solely determined by the initial value F(0). 

The physics governing the escape of F and a from the periodic trap has 
not been discovered. If the value of the energy density g(t) decreases to the 
point where 

lim -A = 

F(t)^0 F(t) 

then lim F ( t )_j. F(t) = 0, and F — > 0; a is controlled by F (because a 3 cosh (F) = 
constant) and will escape the periodic trap and very rapidly transition to a 
constant a(i end f inflation) « oi, (F(t) = Wt > i cnd of in fl at ion)- However 
the governing physics may have more similarity with quantum mechanical tun- 
neling through a barrier than the decay of the controlling parameter g(t) to 
zero. 
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6.1 Dark Photon 



The photon field, carried by a type-2 real bosonic spinor field ip2 a in this 
model, is mapped to the electromagnetic four-vector potential A Q through a 
nonlinear (in u ) interaction with the unit field: 



1 



(TU 



(44) 



A Q is a true SO(3, 1; R) vector field since fl^ is a type-2 spinor. In general, the 
type-2 real bosonic spinor photon field ip2 a is mapped to the electromagnetic 
spacetime potential A# through an interaction with the unit field that includes 
all eight components: 



(Ai, A 2 , A 3 , A 4 ,Ci, 
1 ~ 



-4)D' 



i DB 



D 



aw 



(45) 



The (Ci, C2, C3, C4) are 50(3, 1; R) scalars; for notational convenience, the 
(Ci,C2,C3) are collected into C. 

For the representations employed in this paper 



A 



1 



r 1p2 



e 2 



2 i]j 2 + e 2 ip 2 

-£ , 1 
2 ip 2 

F_ . I F_ . c 

2 ^ 2 — e 2 
2 ^2 - e 2 ^2 

-£ , 6 
2 ^ 2 

6 , 



(46) 



£ / 7 

e 2 ^2 

£ / 2 
e 2 ^ 2 

2 



In a Minkowski limit (F = 0) the electric and magnetic field frame com- 
ponents {E^jB^ 1 } are given (after inflation ends) by 



cdt. 



^A 4 + ^ x C* 



(47) 



and 



^ x^+jfC* + ^C 4 
c -i^C 4 



1 

" 72-g-p 



(48) 
£ = 



c at 

We note that imposing C = and C 4 = — -jt^jC, where ^-^C 
const x magnetic monopolc charge density = removes the (j^,<Cij depen- 
dence from E M and B M . If the (Ci, C2, C3, C 4 ) are decoupled from Standard 
Model degrees of freedom in this manner, then they comprise one set of dark 
degrees of freedom. These dark degrees of freedom receive energy at t = 0, 
just like all of the other degrees of freedom. 
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6.2 Fluctuations of the quantized unit field 

It is widely believed that primordial curvature perturbations sourced by quan- 
tum field fluctuations during the early time history of the universe evolve into 
macroscopic perturbations whose time evolution lead to cosmological structure 
formation. To calculate the quantum fluctuations the bosonic spinor unit field, 
which satisfies nonlinear field equations, must be quantized using an approach 
that guarantees the renormalizability of the quantum theory and the correct 
normalization of the unit field. Quantization should be performed after the g 
and p functions have been redefined and associated to a field, as in Ea.fTB"]. 
Since at comoving time t = a non-vanishing average current drift velocity 
J d exists in this model there will be a non-isotropic term in the total energy- 
momentum tensor. One needs to perturb the FLRW metric for g a p simply to 
obtain the classical solution, then carry out quantization. The investigation is 
ongoing. 



7 Appendix 1: Problem Formulation 

Let X3.1 denote a connected orientable four-dimensional, Hausdorff, smooth 
pseudo-Riemannian manifold whose local tangent spaces are isomorphic to 
flat Minkowski spacetime M 3i i. We assume that the Stiefel - Whitney class 
W2 (X 3j i) vanishes, so that X 3j i admits a spin structure [TT] [TH] and spinor 
fields may defined on this manifold. Let g denote the pseudo-Riemannian met- 
ric tensor on X 3i i. The signature of the metric g is assumed to be (+ + H — ) 
almost everywhere [because the scale factor a(t) may periodically vanish]. 
(X 3j i;(;) is called a spacetime. We assume that g <H> g a p = g a p{x^) car- 
ries the Newton-Einstein gravitational degrees of freedom. Here Greek in- 
dices run from I to 4, and £ M arc local coordinates on X 3j i. We are using 
the usual component notation in local charts. The tensor covariant derivative 
with respect to the symmetric connection associated with the metric g will 
be denoted by a semicolon. A partial derivative is sometimes denoted by a 
comma. When enclosed in parenthesis Greek indices denote flat Minkowski 
SO(3, 1) frame indices; otherwise Greek indices are X 3j i vector indices. Let 
x a = (x, y, z 7 ct) be local comoving coordinates for a neighborhood of p G X 3 ,i. 
Following [5] [TB] we introduce an orthonormal tetrad E^ (with inverse 
E"^) that defines a set of four linearly independent vector fields such that 

E = det (e^J = \/—g 7^ 0. The components arc assumed to be orthonor- 
mal: g a p E?s E^. = t](jj,)(v) = diag(+f , +1, +1, —I) with respect to the X 3j i 
metric g a p. 

In this paper we study a universe that is modeled in terms of a trivial 
fiber bundle X 3j i x R 4 ' 4 with base space X 3j i and fiber R 4,4 , with the trivial 
projection onto the first factor. R 4,4 denotes a real, eight dimensional pseudo- 
Euclidean vector space over R. Let q G R 4 ' 4 and q A G R, A, B, ... = I, 2, 8, 
denote the Cartesian coordinates of q with respect to a definite, right-handed 
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orthonormal frame. 0(4, 4; R) (respectively, SO(4, 4; R)) may be defined as the 
group of all real matrices (respectively, with unit determinant) that preserve 
the quadratic form 

(<Z 8 ) 2 + (q 1 ) 2 + (q 2 ) 2 + (q 3 ) 2 - + {q"f + (q 6 f + (q 7 ) 2 ] ■ 

0(4, 4; R) is a pseudo-orthogonal Lie group that possess two connected com- 
ponents [3] , [T2] , with 50(4, 4; R) being the identity component (the connected 
component containing the identity matrix). Spin(4, 4, R), alternatively denoted 
50(4, 4;R), is the 2-to-l covering gro up of 50(4, 4; R). R 4 ' 4 may be endowed 
with both 50(4, 4; R)-invariant and 50(4, 4; R)-invariant pseudo-Euclidean 
metrics. 

The 50(4, 4; R)-invariant pseudo-Euclidean metric tensor G (respectively, 
inverse G _1 ) may be defined in terms of its components Gab (respectively, 
(G _1 ) = <G AB ) that, with respect to a definite Cartesian, right-handed 
orthonormal frame, are numerically equal to 

GxB = G « = (V._oJ (49) 

where 773.1 = diag(l, 1, 1, -1) is the pseudo-Euclidean metric on flat four- 
dimensional Minkowski space-time M3 1. The indefinite inner product is real- 
ized as R 4 < 4 x R 4 ' 4 9 (a, b) i-> < a ,b >= & A Ba A b B e R. Here A, B, ... = f, 
... , 8 may be regarded as R 4 ' 4 vector indices. 

The fiber R 4 ' 4 also admits a 50(4, 4; R)-invariant pseudo-Euclidean 8x8 
spinor metric a, 



1 

1 



(50) 



Given spinors ip,x € R 4 ' 4 , ipax defines a 50(4, 4; R)-invariant indefinite 
inner product on R 4 ' 4 . Here, following Dirac, ip denotes the transpose of i[>- 

Note that a 2 is equal to the unit matrix, so that the eigenvalues of a of 
are ±1. Since the trace of a is zero, these eigenvalues occur with equal multi- 
plicity. The matrix elements of a are denoted as er a b, a, &, ... = 1, ...8. Gener- 
ally, the early lower case Latin character indices a,b, ... denote Spin(4, 4, R) = 
50(4, 4; R) spinor indices, and run from 1 to 8. A spinor element ip <E R 4 ' 4 is 
defined to have components ip a <G R. The "left-handed" and "right-handed" 
Spin(8) spinors have 50(4, 4; R) counterparts that are denoted ipi and -02 in 
this paper, and transform, respectively, under two inequivalent real 8 x 8 ir- 
reducible spinor representations of 50(4, 4; R) that we call D(x) (type 1) and 
D(2) (type 2). For simplicity, we do not distinguish the spinor index on ip2 from 
that on ipi using a convention such as ipia and ip2 a for spinor components. As 
long as no confusion arises, this simplified notation may be justified. 

When applying this formalism to the physical world we usually restrict the 
allowed transformations as follows: 
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7.0.1 Restriction o/50(4,4;R) to 50(3, 1;R) 

Let (x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ) G R 4 ' 4 . In this paper we assume that for this 
model to possess physical significance the action of 50(4, 4; R) must be re- 
stricted to 50(3, 1;R) in a manner that links transformations of 
to 4 . Therefore we are led to identify R 4 ' 4 with the single particle 

relativistic phase space © *Msa, R 4,4 = M3 i © *Mz,i, and to restrict 
50(4, 4; R) to the subgroup 

K = {<? G 50(4, 471) | g = ( k Q ~° 1 j , lis 50(3, 1; R)| . (51) 

Under the action of L G S0(3, 1; R), x Q i-> x Q = J L Q (3 x /3 and x 4+Q = p a ^ 
Pa = P/9 -^a ! Pa = -Pa a;™ is an invariant, since for g G IK 



7.1 Preferred orthonormal frame 3^ on R 4,4 ; split octonion algebra 

Let Os denote the split octonion algebra. The new "unit field" interaction is 
carried by a distinguished split octonion u G O5 that is algebraically distin- 
guished by the fact that u is the Os multiplicative identity element j2Tj . 
The components of the unit field u with respect to a (non-associative) canon- 
ical oriented spinor basis e a of the split octonion algebra are denoted as u a . 
Throughout this paper we work with the components u a of u, rather than 
with u, to avoid the complexities of the non-associativity of the split octonion 
algebra. In this paper we seek a solution u a to the coupled field equations of 
the form u a = -j= (o, e"^ , 0, 0, 0, , 0, oV 

We do not introduce a characteristic mass for the u field quanta, and 
hence we regard the u unit field to be dimensionless. The real unit spinor 
field u a — u a (x^) carries the massless "unit field" interaction and defines a 
preferred orthonormal frame "$ A on R 4,4 . ^ a A is defined in terms of u and 
natural extension of the Dirac matrices as follows. Define the tau matrices by 
requiring that the tau matrices satisfy (the tilde denotes transpose) 



and 



(TT A = ot a = r A a (53) 

^A—B , ^B—A of irAB -=AB , —B A i KA \ 
T T +T T = 2l 8x8 U ~T T +T T , (54) 



where Isx8 denotes the 8x8 unit matrix. Denoting the matrix elements of r 

a 
b 



A 



by t A a b , we may write Eq. [53] as 



(55) 
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where we have used a to lower the spinor indices. The tau matrices verify an 
important identity [TO] : Let M be any 8x8 matrix satisfying 

oM = oM (56) 



and moreover transforming under 50(4, 4; R) according to 

MhjD^MDf!) -1 (57) 

Then 

t a Mt a =l 8xS tr(M) (58) 
(18x8 denotes the 8x8 unit matrix). 

7.1.1 Components of $ A 

The components $ A of the R 4 ' 4 orthogonal frame are 

F A = —t=T A \u h (59) 
Vucru 

and its inverse is 

^ = —^U C <7 cb T Ab a (60) 

Vucru 

Using Eg. [58] one easily verifies that 

{W8 = «? = «lJi, (si) 

and since a matrix commutes with its inverse, 

si = $2r B . (62) 

Why is this important? The orthonormal frame $ A may be employed to 
map a vector field associated to a fundamental Boson "particle" to a spinor 
field, and then back again. 

7.1.2 Split octonion algebra 

A nonassociative alternative multiplication of the oriented spinor basis e a (re- 
spectively, oriented vector basis e A ) may be defined that endows the real vector 
space R 4 ' 4 with the structure of a normed algebra with multiplicative unit. This 
is accomplished by specifying the multiplication constants m c ab (respectively, 
m AB ) of the algebra, which verify 

6a 66 = e c fn c ab 

e A e B = e c m c AB (63) 
The set of m c ab (respectively, mF AB ) defined by 

m ab = Sa T A ^ 

m C AB =^r A \t B (64) 
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a 




Fig. 1 Dynkin diagrarh for D4 

is a realization of the multiplication constants of the algebra. It has been 
proven that the nonassociative product defined by Eq.[63] and Eq.[64] of the 
spinor basis e a (respectively, of the vector basis ca ) endows the real vector 
space R 4,4 with the structure of the split octonion algebra over the reals [20j . 
The multiplicative unit element of this split octonion algebra is 



8 Appendix 2: Spinor representations of SO(4, 4;R) on R 4 ' 4 
8.1 Representations of SO(8,C) 

Both 0(4, 4; R) and 5*0(4, 4; R) are real forms of the classical complex or- 
thogonal group 0(8, C). The relationship between Clifford algebras C„ and 
the spinor representations of the classical complex orthogonal groups is well 
known. Summarizing a few basic facts that may be found in, for example, 
Boerner, The Representations of Groups [3], the Clifford algebra C% may be 
defined as the algebra generated by a set of 8 elements ej,j = 1, . . . , 8, that an- 
ticommute with each other and have unit square &j + e^ ej — 25 jk Ii6x 16, I 
= unit matrix. The scaled commutators j (ej — Sj) computed from the 
16-dimcnsional irreducible representation of the ej are the infinitesimal gen- 
erators of a reducible 16-dimcnsional representation of Spin(8,C), which is 
the universal double covering of the special orthogonal group S0(8,C). This 
16-dimcnsional representation of is fully reducible to the direct sum of two in- 
equivalent irreducible 8x8 spin representations of the infinitesimal generators 
of Spin(8,C), which leads to the identification of type 1 and type 2 spinors. 
The fundamental irreducible vector representation of S0(8,C) is also 8x8. 
The Dynkin diagram for D4 = SO (8) is symmetrical and pictured in Figure 
1. The three outer nodes correspond to the vector representation (left- most 
node), type 1 spinor and type 2 spinor representations of Spin(8), and the 
central node corresponds to the adjoint representation. 

As is the case with S0(2n;C), the basic spinor representation of the 
pseudo-orthogonal group SO(A, 4; R) may be constructed from the irreducible 
generators t A ,A = 1, 8, of the pseudo-Clifford algebra C44 [3], [4], [15]. Lord 
[15] has indicated a general procedure for constructing the spinor representa- 
tions of SO(2n;C) of the first and second kind (handedness) from the genera- 
tors of C2n-2- Following Lord we shall call such irreducible C2 n -2 generators 
"reduced Brauer-Weyl generators" [3] . Using Lord's procedure we construct 
the representation of 50(4, 4; R) by first defining the real 8x8 matrix reduced 



id = e a u a = es = u. 



(65) 
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Brauer-Weyl generators t a ,t a 7 A,B,... = 1,...,8, of the pscudo- Clifford al- 
gebra that anticommutc and have square ±1 according to Eq.[S3] and 
Eq. [54] . We adopt a real irreducible 8x8 matrix representation of the tau ma- 
trices (see the Appendix) in which r 8 = Igxs = T§ - Then by Eq.[M] T A = —t a 
for A = 1, ... ,7. Hence, again by Eq.[S3], ir A ) 2 is equal to — I 8X 8 f° r A = 1,2,3 
and is equal to +Isx8 f° r A = 4,5,6, 7, 8. 



8.2 Transformation under action of £0(4, 4; R) 

The special Lorentz transformation properties of the theory may be deter- 



mined by constructing a real reducible 16 x 16 matrix representation of £0(4, 4; 1 
utilizing the irreducible generators t , A = 1,...,8 of the (pseudo-) Clifford 
algebra C4.4. Continuing to employ Lord's general procedure [TS] we define 
the irreducible generators t A as 



Let g G £0(4, 4; R). The 16 x 16 basic spinor representation of £0(4, 4; R) is 
reducible into the two real 8x8 inequivalent irreducible spinor representations 
D(\){g) and D^(g) of £0(4, 4; R). The reduced generators of the two real 
8x8 spinor representations D^(g) and D^(g) of £0(4, 4; R) follow from the 
calculation of the infinitesimal generators 

't a t b -t b t a 

T A T B - T B T A 



t A t B _ t B t A = 



D; • AB 



(1) 

"(2) 



D r AB 



(67) 



of the 16-componcnt spinor representation of £0(4, 4; R). We see that, as is in 
fact well known from the general theory, the 16-componcnt spinor representa- 
tion of £0(4, 4; R) is the direct sum of two (inequivalent) real 8x8 irreducible 
spinor representations Dn\ = D^(g) and D(2) = D/ 2 \(g) of £0(4, 4; R) 3 g 
1 {1) AB and D {2 



that are generated by Dii\ AB and Dt 2 ) AB respectively. This we record as 



4D {1) AB =t a t b -t b t a (68) 

and 

4D {2) AB =t a t b -t b t a (69) 

For completeness we remark that the generators of the two spinor types 
are images of the projection operators 

X ± =4(l±i 9 ) 



x+ = 



2 

I 8 X8 




0I 8 x8 



(70) 



21 



Patrick L. Nash 



where 

t 9 = tHHHHHHh 8 = ^ t ° ) . (71) 

Here 

=0 — 12— 3 4— 5,fi— 7 1 2 3 4H fi^7 / 7 o\ 
T =TTTTTTT~TTTTTTT {> ^) 

and 

T =TTTTTTT = — T T T T T T T = — T [ I O ) 

The representation of the tan matrices is irreducible, r has square equal 
to +Isx8 an d commutes with each of the t a matrices (and therefore with all 
of their products). Therefore we conclude that r° = ills x 8 m any irreducible 
representation. r° = I§ X 8 m the irreducible representation given in the Ap- 
pendix. 

Let ujab = —uba £ 8, 4, 5 = 1, . . . , 8, enumerate a set of 28 real parame- 



ters that coordinatize g = g(uj) e 50(4, 4; R). Also, let L = L(g) <G 5<3(4,4;R) 
have matrix elements L A B , denote the real 8x8 matrix with matrix ele- 
ments lu a b = Gr Ac u>cB, ui = \loabD(i) AB and lo 2 = \ujabD( 2 ) AB . We find 
that 

D(i) = D {1) (g) = cxp 
D( 2 ) = D(2){g) = cxp 

L A B = L A B (g) = {cxp(J)} A B (74) 



where, under the action of 50(4, 4; R), 

Df 1} AB <J = -aD {1) AB => D {l) a = aD^- 1 (75) 
d{ 2) AB <j = -aD {2) AB => D (2) a = aDpf 1 (76) 



J C^AB-Li D — D — 1 J^^-Lj f C D 



B^ d 


— Grc D 1 


= |lgl 


A T B 
B T 




r A D (2) 


A T B 

' B T 




r A D (1) 



(77) 

L A B r B = D {1) -'r A D {2) (78) 
L A B r B =D [2) - 1 r A D [l) (79) 



The canonical 2-1 homomorphism £0(4, 4; R) — > 50(4, 4; R) : g i-» L(g) is given 

by 

8 L A B = tr (D^- 1 r A D {2) t c ) G CB , (80) 

where tr denotes the trace. Note that D^(g(cu)) = D^(g(u>)) when u>as = 0, 
i.e., when one restricts 50(4, 4; R) to 

50(3, 4; R) = {g E 50(4, 4; R) | 

exp (\ujabD{i) AB ) \ \ . . 

n (l n AB\ ) andui AS = 0> (81) 

exp ^ j4B f(2) )/ J 
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[one of the real forms of Spin(7, C)]. When stating mathematical spinor rela- 
tionships we almost always restrict 50(4, 4; R) to 5*0(3, 4; R) without explicit 
notice. 

We adopt a real irreducible 8x8 matrix representation of the tau matrices 
(see Appendix ITU]) in which r 8 = I 8x8 = t 8 . Then by Eq.[33] f A = -t a for A 
= 1, ... ,7. Hence, again by Eq.[M], {t a ) 2 is equal to — Isx8 f° r A = 1,2,3 and 
is equal to +1.8x8 for A = 4,5,6, 7, 8. 



8.3 50(4, 4; R) covariant multiplications 

Let Vi,V2, and V3 be vector spaces over R. A duality is a nondegencrate bilinear 
map V\ X V2 — > R. A triality is a nondegencrate trilinear map V\ x V2 x V3 — >• R. 
A triality may be associated with a bilinear map that some authors call a 
"multiplication" [2] by dualizing, V\ x V2 — >• * V3 = V3. There are two especially 
simple multiplications that possess covariant transformation laws under the 
action of 5*0(4, 4; R). Assume that u i-> u = D^) u, ^1 4 ^ = Dm "01 an d 
V"2 !->• V J 2 = ^(2) "02 under the action of 50(4, 4; R). The first multiplication 
mi A : R 4 > 4 x R 4 ' 4 -> R 4 < 4 is defined by 

Q A = u a t a ip 2 . (82) 

For fixed u a , Q A € R 4 ' 4 depends on 8 real parameters arranged into the type-2 
spinor ip%. 

The second multiplication m 2 AB : R 4 ' 4 x R 4 ' 4 -> F 3 has an image m 
V3 = R 4 ' 4 x R 4 ' 4 , and depends on 8 real parameters (for fixed u a ) arranged 
into the type-1 spinor tpi: 

Q AB =u a r A t b ip\. (83) 

We emphasize that, for fixed u, Q AB possesses only 8 degrees of freedom 
corresponding to the 8 independent degrees of freedom of tpi , so we also refer 
to this map as a "multiplication." Eq.[53] may be solved for the components 

V>i a = i>i a (Q AB )- 

Suppose we break the 50(4, 4; R) symmetry down to 50(3, 1; R) according 
to paragraph 17.0.11 We note that G is invariant under 50(3, 1;R) since it is 
invariant under 50(4, 4; R). It is convenient to define a 50(3, 1; R)-invariant 
symplectic structure Q on R 4 ' 4 (and a complex structure on the split octonion 
algebra) by 

»-(-°.S) 

where denotes the 4x4 zero matrix and 1 denotes the 4x4 unit matrix. 

qAB ma y b e expressed in terms of a linear combination of 'simple' objects 
that transform as M% t i tensors of appropriate rank. Consider first two arbitrary 
-^3,1 vectors with components A M and B M . Q AB may be written as a linear 
combination of A Q M t3 ±A 13 M a , e Q ^A M B„ and rf 3 r/^A^ M v , but this is does 
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not yield a linear relationship between Q and 



t }. Also, since Q 



AB 



has exactly eight independent degrees of freedom it cannot be expressed in 
terms of a traceless symmetric rank 2 tensor, which has nine independent 
components. However we find that the Q AB may be represented in terms of an 
arbitrary antisymmetric M 3y i rank 2 tensor F' 3 " = — ¥ a/3 and two 50(3, 1;K) 
scalars q^ = B4 and q% = £4 according to 



n AB_ I 



AS 



AB 



(85) 



. 1 f a0fj,va 



a/3 



. Note that 



where *F Q @ is dual to F a/3 and defined by 
q[.ab] _ 1 (q^ b _ Q BA ^ is independent of q 8 , that q$ may be expressed as 
Q 88 , and that q± may be expressed as Q 48 . 

If we employ the representation of the tau matrices given in the Appendix, 
and then solve for the Vi a for the case u = ^=(0,1,0,0,0,1,0,0) then we 
obtain 

/-B3 + BA 



^ 1 = 7f 



B x - 
-B 3 



E4 
E 2 
Ei 
-B 4 
-Ei 
E 2 



(86) 



\B2-Bh. J 



9 Appendix 3: Spin Connection 

Recall that Greek indices run from 1 to 4. When enclosed in parenthesis they 
are flat Minkowski spacetime SO(3, 1) indices; otherwise they are ^3,1 vector 
indices. Thee metric tensor g a p has signature (+ + H — ). Let x a = (x,y,z,t) 
be local comoving coordinates for a neighborhood of p G X3 t i. 

Following [5] [TB] we introduce an orthonormal tetrad Ea (with inverse 
E-Qjj) that defines a set of four linearly independent vector fields such that 

E = det (e^?^ = \f—g 7^ 0. The components are assumed to be orthonormal: 

9afi Efy E^ v) = J7( A ,)(„) with respect to the X 3) i metric g aj3 . 
Let 

7 Q = % >V 7 (87) 

The Lorentz connection 0$*'^ and spin connection = gUj,( a )(m py( a )"y(0) —jWj( a )^ 
are defined [TB] using 

-2 W W(« = - - S(« K - d K E^) 

+ E^*E^° (d a E™ - d K E^) E Wfl . (88) 
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Note that the covariant derivative V^g 7™ =7 Q ( g = dpl a + {^",9} r ' y +7° ~ 

7° = 0. Here i s the Christoffel symbol of the second kind. For brevity 

we write 7" V Q u = 7° [u. Q — r a u] = 7" it| Q . 



10 Appendix 4: Irreducible Representation of the r Matrices 

The particular irreducible representation of the tau matrices employed in the 
calculations is 
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(89) 
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-4-2 2 



2(Vcosh(F,-) +l) 

1*1 £ 

2 V cosh(F,) 
Singularities at t = ± 

2(Vcosh(F,) ±l) 
Possible initial states at t = ± 

Fig. 2 -F(i) (blue) Fj (red, dashed) [limiting case of no Fermions and zero electromagnetic 
fields]. [Color online] 
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